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TECHNICAL ROTE 5779 

TENTATIVE MESEHOD FOR CALCDLATIOR OF THE SOUND FIELD ABOUT 
A SOURCE OVER GROUND CONSIDERING DIFFRACTION 
AND SCATTERING INTO SHADOW ZONES 


David C. Pridmore-Bravm and Uno Ingard 


SUMMARY 


A semlenglrlcal method is given for the calculation of the sound 
field about a source over ground considering the effects of vertical 
temperature and wind gradients as well as scattering of sound hy turhu- 
lence into shadow zones. The diffracted field in a wind-created shadow 
zone is analyzed theoretically in the two-dimensional case a n d is sliown 
to he similar to the results obtained for a temperature -created shadow 
field as given in NACA TN 3k9h, The frequency and wlnd-velocilgr depend- 
ence of the scattered field into the shadow zone is estimated from 
LighthlU's theory, and on the basis of these two field contributions 
a semien5)irical formula is constructed for the total field which con- 
tains two adjustable i>arameters. From this expression a set of charts 
has been prepared showing equal sound-pressure contours at 10 feet above 
ground for various source heists, wind velocities, and frequencies. 

The two adjustable parameters in the formula were obtained from 
measTirements using a relatively small source height (lO feet). The 
parameters should actually be a function of heigpit determined by the 
wind Eind temperature profiles. However, in these preliminary calcula- 
tions the parameters have been kept constant, and the fields, particu- 
larly for large source heists, must be considered as preliminary esti- 
mates to be corrected when more information is available. 


IHTR033UCTI0N 


As is well known, sound energy follows curved paths in the presence 
of a temperature gradient. Theory indicates that near a sound source 
this phenomenon has only a sli^t effect on the intensity field, the 
correction being of the second oidLer in the temperature gradient. How- 
ever, if a sound source is placed above ground in a negative tenperature 
gradient (lapse rate), a shadow zone will form at a certain distance 
from it because the rays are screened off by the ground from entering 
this region. The same phenomenon occurs in the presence of a wind 
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gradient^ a shadow heing formed on the tgnwind side of a source If the 
wind velocity increases with hei^t. 

Previous reports on the present pro;ject in atmospheric acoustics 
have described two sets of measur^nents, one performed in the laboratory 
(ref. 1) and the other in the field (results as yet unpublished). In 
the laboratory measurements a uniform (static) tenperature gradient was 
set up in a tro-dimensional propagation chamber. Below a certain fre- 
quency this chamber allowed cylindrical waves to propa^te between two 
parallel steel plates spaced a distance apart less thaq ba-lf a wavelength. 
The temperature gradient was hi^ enou^ to secure a sound shadow within 
a short distance from a source placed near one end. The measured shadow 
field in the chamber agreed well with diffraction theory, which predicts 
that the sound pressure in this region should fall off eacponentiaUy with 
distance from the shadow boundary at a rate which increases with the fre- 
quency and the tenperature gradient. Finally, the two-dimensional analy- 
sis of the diffracted field in a wind-created shadow, which is included 
in the appendix, predicts the same behavior of the sound field as was 
predicted for the temperature-created shadow. 

Meastirements in the field taken in a wind-created shadow gp.ve 
results of the type illustrated qualitatively in figure 1. The sound 
intensity at a point within the shadow was observed first to diminish 
with frequency in accordance with diffraction theory, but after a cer- 
tain frequency it increased again, thus indicating the presence of 
another mode of energy penetration into the shadow which predominated 
over diffraction at the hl^er frequencies. It seems reasonable to 
stppose that this contribution to the hi^er frequencies is due to the 
presence of turbulence in the atmosphere which scatters sound from the 
normal zone into the shadow as illustrated schematically in figure 2. 

On the basis of this mechanism a semienpirical theory is given 
which attempts to eaqplain the field measurements and to produce an 
approximate method by which the sound field can be predicted. 

The present investigation was conducted at the Acoustics Laboratory 
of the Massachusetts Institute of Technology under the sponsorship and 
with the financial assistance of the National Advisory Committee for 
Aeronautics. 


SYMBOLS 


A,B enplrlcal constants 

c speed of sound, fps 

frequency, . cps 


f 
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H relative Irumidlty, percent 

h. sotirce hei^t, ft 

I sound intensity 

k propagation constant, 2rtf/c 

p sound pressui« field 

R source-to-receiver distance, ft 

Rq distance between source at hei^t h and receiver at 

height 10 ft, ft 

r radial coordinate, distance from source, ft 

rQ horizontal component of Rq, ft 

s wind-velocity gradient, ^ ^ 

t time, sec 

V wind velocity, fps unless otherwise indicated 

X horizontal distance from source, ft 

Xq horizontal distance from source to shadow boundary, ft 

z vertical coordinate, receiver hei^t, ft 

a propagation factor 

P = v/c 

T temperatiare gradient, °C/ft 

angle between wind and sound direction^ veloci-ty potential 

0 ) = 2jrf 

AMIYSIS 
Diffracted Field 

A two-dimensional calculation of the sound field from a source 
located in a wind gradient, given in the appendix, predicts that within 
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the acoustic shadow on the upwind side of the source the pressure field 
due to diffraction alone should he of the form 

p ~ exp I^B * (x - 

where Xq is the horizontal distance to the shadow-zone hovtndary, 

k = 2itf/c, s = (l/c)(dv/dz), and B* = O .7 for a hard houndaiy. In 
three dim^sions it seems probable that the same relation will hold in 
the direction against the wind with the addition of the factor l/V? 
to take account of divergence. Moreover, in practice it turns out that 
the wind velocity v and its gradient are roughly proportional. Sub- 
stituting one for the other and at the same time replacing B* by B 
to allow for the unknown proportionali-ty leads to the assumption for the 
diffracted sound in the shadow zone of an expression of the form 

Ig ^ g-2B(r-ro)fl/5v2/3 


where Ip is the sound intensity on the shadow boundary. It is 
assumed that Iq varies Inversely as the square of the distance from 
the source, that is, that it has the form where Eq is the dis- 

tance from the source at hei^t h to the shadow boundary at hel^t z. 


Scattered Part 

In addition to the sound field diffracted into the shedow region 
there will also be a contribution from the scattering of sound from 
turbulence in the sound region, as already indicated schematically in 
figure 2. The theory of the scattering of soTxnd from turbulence pre- 
sented by Li^thill (ref. 2) is essentially applicable to only a small 
region of turbulence in space and has to be modified to a multiple - 
scattering theory in order to apply in the present case. However, with- 
out having performed such a detailed analysis, it is assumed that in 
this case the frequency and velocily dependence of the scattered energy 
will still be the same as that given by Lighthill. In particular, atten- 
tion is limited to that region in which the sound wavelength can be con- 
sidered small conpared with the scale of the turbulence. In that case 
the scattered energy should be proportional to the sound frequency squared 
as well as to the mean square of the velocity fluct\xations. Since the 
velocity fluctuations are proportional to the average wind velocity v, 
the scattered sound may therefore be set proportional to f^v^. Although 
the scattered sound must originate from all over the sound zone, in this 
simple analysis it is assumed that this distributed source can be replaced 
by a single source located between the receiver and the shadow boundary. 
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It is also assumed that the sound intensity from this equivalent source 
spreads out in a spherical wave and accordingly may he written as pro- 
portional to where R is the source-to-receiver disteince. 

To insure that this equivalent sound source lies within the normal, zone 
a is set equal to ^ horizontal distance to the 

shadow boundary. The exact value assigned to a turns out to be unim- 
portant since, for a large value of R, it will have a negligible effect 
on the VEilue of the term, and for a small value of R this term is 
itself small compared with the one representing the diffracted sound. 


Total Field 

From the previous discussions the following expression is obtained 
for the total sound intensity within the acoustic shadow: 


I = 



g-2B (r-ro) f l/5v2/3 


( 1 ) 


The assumption is made here that the shadow distance ro is related to 
the source heiggh t h as it is for a linear wind gradient; namely, 

Tq - Kv"^/^(l + l/h/z j where K is a constant. The receiver hei^t z 
is put equal to 10 feet throughout the calculations. From the defini- 
tions it follows also that R^ = r^ + (z - h)^ and Rq^ = + (z - h)^. 

The unknown constants A and E can be found by fitting equation (l) 
to data obtained in field measvurements. 

The constant A is essentially a measure of the relative strengths 
of the diffracted and scattered sound pressure levels and for the pur- 
pose of this analysis is considered to be ind^endent of the other phys- 
ical quantities. Similarly, B is assumed to d^end only on the ground 
impedance. The constants A and B can be estimated empirically from 
a knowledge of (l) the critical frequency at which the souctd pressure 
level at a given point in the shadow is a minlrmnn under a known set of 
physical conditions an d (2) the total attenuation brought about Ty wind 
effects, which is measured under these same conditions. Row, in one 
set of field measxirements (ref. 5) made under the following conditions 

r = 1,000 ft upwind of source 

ro = 200 ft 

V *= 15 fps 


h = 10 ft 
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the critical frequency turned out to he about 1,000 cps. At this fre- 
quency a total attenuation of 37 decibels more than the val\ae at the 
shadow boundary was measured; in other words, from equation (l). 


10 log 


(afV/ro )2 

+ (a 



[fy(r - J 4 

2 

+ 

(A/rro) 

4 


57 


( 2 ) 


where 7 = B(r - rQ )fl/V/3. 

Differentiation of equation (l) with respect to frequency gives A 
in terms of 7 and critical frequency as 



Substituting for A in equation (2) gives 7 = l^.l nepers whence 

A s 50 ft"^ sec“^ 

B = 8 X 10“5 sec ft”5/5 

Finally, the relation adopted between shadow distance and wind veloc- 
i-ty Tq = Kv"^/^(l + ^h/zj determines K as 

K = 400 ft5/2 sec“V2 
KOMERICAL CALCULATI0H5 


Equation (2) giving the critical frequency is plotted in figure 3 
against wind velocity for several values of the source-to-receiver 
parameter r and for a so-urce hei^t of 10 feet. These curves divide 
the upper ri^t-hand region where scattering predominates from the 
lower left-hand region where diffraction is inportant. The so\md 
intensity given in equation (l) in terms of distance, frequency, wind 
velocity, and source hei^t has been plotted as a function of distance 
for several values of the other parameters. Figures 4(a) to 4(f) are 
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plots of the soilnd-pres sure-level Increment (or 10 log (I/Iq) where 
lo = Constant) on the igjwlnd side of a source located 10, 100, and 
1,000 feet above ground and driven at 2^0 and 500 cps in various winds. 
From these plots polar diagrams giving the equal sound-level contours 
aroxind the source were prepared iinder the following two assumptions: 

(l) Within the shadow, formula (l) was assumed to hold for an angle 0 

between wind and sound provided that v was r^laced by v cos 0,^ 
and (2) within the normal region, the sound intensity was assumed to fall 
off with the square of the distance from the soiirce. These polar diagrams 
are given in figures 5(a) to 5(2) for wind velocities of 10 and 20 miph, 
frequencies of 250 and 500 cps, and source heists of 10, 100, and 
1,000 feet wiftic-fTig a total of 12 diagrams. In each of these the position 
of the shadow boundary has been drawn in as a dashed line. The decibel 
values labeling the equal sound-level contours Indicate the decrease of 
sound pressure level below the reference sound pressure level as obtained 
in a free field at 100 feet from the source. 

As a specific example of the use of these diagrams, consider a sound 
source such as a jet to be located 10 feet above ground in a 20-mph wind 
measured at 10 feet above ground. The problem is to predict the sound 
pressure level 10 feet above ground at a distemce of 1,000 feet from the 
source. Assume the power output of the engine to be such that at 100 feet 
its sound press-ure level is 120 decibels and assume that it is mainly con- 
centrated around 250 cps. The numbers given in the charts Indicate the 
mmiber of decibels that has to be subtracted from the value at 100 feet 
in order to obtain the sound pressure ' level at the point in question. 
Therefore, in this example the sotind pressure level on the downwind side 
should be approximately 120 - 20 = 100 decibels at 1,000 feet, while on 
the upwind side it becomes 120 - 46 = 7^ decibels at that distance. In 
the direction 4-5° into the wind the level should be 80 decibels. 

The formula presented in equation (l) can, of course, be tised to 
prepare charts for more extended reuages of the variables ^considered. 

If, in addition to the wind gradient, there also exists a temperature 
grsuiient t, the quantliy v in equation (l) must be replaced by 
V cos 0 + 10^ (t is the temperature gradient in °C per foot). Thus, 

if V < 10^, at some angle 0 given by cos 0 = -v cos 0/loV the 
effects of wind and temperature gradients cancel each other.- This angle 
demarcates the boundary between the shadow and the normal zone. The 
charts presented herein refer to the situation when there is no tempera- 
ture gradient (t = 0°/ft and 0 = 90°). 


^A preliminary analysis of the three-dimensional wind equation 
indicates that this procedure is justified. 
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It should he noted that from the sound pressure level obtained here 
should he subtracted the effect of air absorption which can he written 
as approximately 



- ■ ■ db/mile 

20 + H 


where H is the relative humidity in percent. This effect is negligible 
in most cases. For exampj.e^ at 250 cps and 50-percent relative humidiiy 
it amounts to an attenuation of only 0.8 decibel in a mile while at 
500 cps the attenuation is 2.3 decibels per mile. 


Massachusetts Institute of Technology, 
Cambridge, Mass., A^ugust 19 , 1955 * 
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APPENDIX 

TWO-DIMENSIOML AMLYSIS OP DIFFRACTED FIELD IN 
WIND-CREATED SHADOW 
Cylindrical. Source in Wind Gradient 


Consider a cylindrical, source to be located a hei^t b above 
ground in a horizontal wind. The cylinder is parallel to the ground 
and perpendicular to the wind direction. The sound field at large dis- 
tances from such a source is sou^t as the solution to the wave equation 
in a moving mediomi (ref. which is 


= 



0 


If the x-axis is taken in the wiM direction and a time factor e“ 
is canceled the equation becomes approximately, for p(z) = v/c « 1, 


(v2 + + 21k^ » 0 


Two-dimensional solutions to this equation are of the form 

P(z)eC“ 

where 

^ + (a2 + k2 + 21kpa)F = 0 

dz^ 

For simplicity a linear wind gradient ^ = -sz is assumed (s constant) 
for which 


F(z) . (v) 
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and is tlie one-third order Eankel function of the first kind. 

If these solutions are reqijlred to satisfy the houndary condition 
dP/dz = 0 at z =s 0 (coi^esponding to infinite ground 'liqpedance), an 
infinite set of functions Fji(z) with corresponding propagation fac- 
tors Oj;^ is obtained. The sound field nmy he e^cpressed as a sum of 
these functions in the form 


0 = 21 1‘nOi)Fn(z)e‘^ (A2) 

n 

The houndary condition dP/dz = 0 at z = 0 leads to 

= 0 

whence, 

kos 


where B = 0.685 for the first root. The propagation factor a is 
then 


a, = tikjl - (3sk:"2aB)^'^^e^’^/^ 

Here, in order to insure outgoing waves in conjunction with the 

time factor e"-^^, one must vise the plus sign to the ri^t of the 
source and the minus sign to the left. Since (for high frequencies) 
the second term in the brackets is small compared with the first, one 
may put approximately ct = ±ik in this tbrm and so obtain for a on 
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expanding (a) to the rl^t of the source. 



and (h) to the left of the source. 



To the left of the source , a is purely Imaginary and leads to 
undan^ed propagation, whil e to the right It Is con^>lex, indicating that 
the waves are attentaated in this direction. This is to he expected 
from the form of the velocity function which leads to an effective 
refraction of the waves, downward to the left and, upward to the ri^t. 


Hl^-Frequency Behavior of Solution 

The asyuptotic form of eaeh mode of .the solution (A2) for hi^ fre- 
quencies is (apart from a phase factor) 


0 ^ ^i Jw(z)iw(h)4<xx-<w^ 


where w(z) is given in equation (Al). If it is assumed that 
a = ik(l - 6), where from equation (A3) B = on the 

upwind side of the source, then to the first order in B and sz 


[k^ _ k2(i _ 26) + 


3k2s • 



Now, from ray theory the horizonteil distance from the source at 
hei^t h to the shadow houndary at height z is 
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while the corresponding travel time of the ray is 

to - # (if ''^(i + f th) 

Rewriting w in terms of these parameters gives 

w(t) + w(h) = k(xo - ct^) - ]£Xq6 , 

The hi^-frequency behavior of the solution is then 

^ ^-i [w (t-to) -h (x-Xo) (l-B^ 

Within the shadow zone (x > Xq^ the solution is essentially represented 
by the first mode, which is attenuated at a rate given by 

k lin 8 = -^(5Bs)^/^k^/5 

**» 0.7 s^/^k^/5 neper s/distance 
where Im signifies the imEiginary part. 
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Figure 1.- Aa illustrative plot of frequency spectrum of sound intensity 
in shadov showing departure from diffraction theory. 


WIND 



SCATTERING OF 
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SOUND ZONE 
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SHADOW ORi 


Figure 2.- Schematic diagram illustrating how sound penetrates shadow- 
owing to the presence of turbulent scattering in sound zone. 
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Figure 3*- Plots of critical frequ^cy as function of wind velocity for 
several values of source-to-receiver parameter r for source 
hei^t of 10 feet. Curves divide upper ri^t-hand region where scat- 
tering predominates from lower left-hand region where diffraction 
predominates . 
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(a) Source hel^t^ ID feet; frequency, 250 cps. 

Figure 4.- Increment of sound paressure level as a function of distance 
from, source for various vlnd velocities. Wind velocity Is measured 
10 feet above ground and Is assumed to increase vltb. hel^t; ordinate 
gives increment of sound pressure level from reference sound pressitre 
level as obtained In free field at 100 feet fron source. 
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t, 100 feet; frequency, 250 c^s 
ure 4 .- Continued. 
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(c) Source height, 1,000 feetj frequency, 250 cps. 
Figure 4.- Continued. 
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(e) Source height, 100 feet} freguencyi 500 cps 
Figure 4 -.- Coirtinued. 
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(f) Source hei^t, 1,000 feet; freq.uency, 500 cps. 
Figure 4.- Concluded. 
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(a) Source hei^t, 10 feet; wind velocity, 10 ' b®1x; frequency, 250 cps. 

FlfTure 5.- Equal sound-level contours measured 10 feet atove 
^out source in wind for different wind velocities, source heights, 
and frequencies. Wind direction is toward 0°; decltel values on 
contours indicate decrease of sound pressure level Taelow ohe ref^- 
ence sound pressure level as obtained in free field at 100 feet from 
source; dashed line Indicates position of shadow boundary. 
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(j) Source hei^t, 100 feet; wind velocity, 20 mph; frequency, 500 cps 


Figure 5 .- Continued. 
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(k) Source hei^t, 1,000 feet; wind velocity, 10 mph; frequency, 500 cps 

Figure 5.- Continxied. 
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(Z) Source height, 1,000 feet; wind velocity, 20 raph; frequency, 500 cps. 

Figure 5«- Concluded. 





